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Electronic and magneto-optical properties of monolayer phosphorene quantum dots 


Rui Zhang\ X.Y. Zhou^ D. Zhang\ W. K. Lou\ F. Zhai^, Kai Chang^* 

^SKLSM, Institute of Semiconductors, Chinese Academy of Sciences, RO. Box 912, 100083 Beijing, China and 
^Department of Physics, Zhejiang Normal University, Jinhua 321004, China. 

We theoretically investigate the electronic and magneto-optical properties of rectangular, hexangular, and 
triangular monolayer phosphorene quantum dots (MPQDs) utilizing the tight-binding method. The electronic 
states, density of states, electronic density distribution, and Laudau levels as well as the optical absorption 
spectrum are calculated numerically. Our calculations show that: (1) edge states appear in the band gap in all 
kinds of MPQDs regardless of their shapes and edge configurations due to the anisotropic electron hopping in 
monolayer phosphorene (MLP). Electrons in any edge state appear only in the armchair direction of the dot 
boundary, which is distinct from that in graphene quantum dots; (2) the magnetic levels of MPQDs exhibit a 
Hofstadter-butterfiy spectrum and approach the Landau levels of MLP as the magnetic field increases . A ’’flat 
band” appears in the magneto-energy spectrum which is totally different from that of MLP; (3) the electronic 
and optical properties can be tuned by the dot size, the types of boundary edges and the external magnetic field. 


PACS numbers: 73.21.-f, 78.67.-n, 75.75.-c, 81.07.-b 


I. INTRODUCTION 

Since its discovery in 2004 \ graphene has become a ma¬ 
jor issue of theoretical and experimental study in condensed- 
matter physics However, graphene suffers from limitations 
for practical applications due to the absence of a band gap. 
The discovery of graphene evokes great interest and triggers 
intensive investigation on two-dimensional (2D) atomic-layer 
systems, such as monolayered hexagonal BN, transition metal 
dichalcogenides (TMDCs)"^"^, silicene^’^, and germanane^’^. 

Recently, a new 2D material, few-layer black phospho¬ 
rus (termed ’phosphorene’), was successfully fabricated from 
bulk black phosphorus in experimentBP consists of 
phosphorus atom layers, where adjacent layers coupled by 
weak Van der Waals interactions^^. Inside a single layer, 
each phosphorus atom is covalently bonded with three nearest 
phosphorus atoms to form a puckered honeycomb structure 
(Fig. 1(a)). BP possesses a direct band gap (0.3 eV) located 
at Z point. This direct gap increases to 1.5-2 eV when the 
film thickness decreases from bulk to few layers and eventu¬ 
ally monolayer via mechanical exfoliation. This band gap is 
crucial for practical applications in field-effect transistors^^. 
Unlike other widely studied 2D materials such as graphene 
and TMDCs, the band structure, electrical conductivity, ther¬ 
mal conductivity, and optical responses of few-layer BP are 
all highly anisotropic^ 

As a newly emerged member of the 2D crystal family, 
phosphorenes have recently gained significant attention due 
to their interesting electrical and optical properties. To date, 
various properties of phosphorene have been predicted theo¬ 
retically and verified experimentally, including those related 
to strain induced gap modification^^, field transistor effect^^, 
anisotropic Landau levels^^, quantum oscillationsand 
quantum Hall effect^^. An important aspect on the applica¬ 
tions of phsophorene is the finite size effect, which has been 
explored mainly in phosphorene nanoribbons^^”^^. Phospho- 
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rene quantum dots have not yet been researched extensively 
so far. 

Motivated by the recently experimental sysnethis of few- 
layer phosphorene quantum dots by chemical method^^, we 
investigate the electronic and magneto-optical properties of 
monolayer phosphorene quantum dots (MPQDs) here. The 
electronic structure of monolayer phosphorene (MLP) can be 
well described by a tight-binding (TB) model^^. Five types of 
MPQDs with a rectangular, hexangular, and triangular shape 
are considered. We find edge states detached from the bulk 
band gap in all kinds of MPQDs studied regardless of their 
shapes and edge configurations due to the anisotropic hop¬ 
pings parameters in MLP. The magnetic levels of MPQDs ex¬ 
hibit a Hofstadter-butterfiy spectrum and approach the Landau 
levels of MLP as the magnetic field increases. A ’’flat band” 
appears in the magneto-energy spectrum which is totally dif¬ 
ferent from that of MLP The electronic and optical properties 



FIG. 1: (Color online) (a) Top and (b) side view of phosphorene. The 
hopping terms in the tight-binding Hamiltonian are shown by edges 
between sites. There are four phosphorus atoms in a unit cell labeled 
A, B, C, and D. 
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can be tuned by the size and edge type of MPQDs and the 
external magnetic field. 

The paper is organized as follows. In Sec. II, we present 
the TB model. In Sec. Ill, we calculate the electronic proper¬ 
ties including the density of state (DOS), probability distribu¬ 
tion of edge states, and size dependent energy gap of MPQDs. 
In Sec. IV, we investigate the magneto-energy spectrum and 
magneto-optical spectrum. Finally, we summarize our results 
in Sec. V. 


MPQDs (HAMPQDs), hexangular zigzag MPQDs (HAM- 
PQDs), triangular armchair MPQDs (TAMPQDs) and trian¬ 
gular zigzag MPQDs (TZMPQDs). The basic parameters of 
MPQDs and the information of five typical MPQDs are shown 
in TABLE. I. The total number of phosphorous atoms in a 
MPQD, Ntotah is determined by the number of hexagonal units 
along the armchair and/or zigzag direction {Na and/or N^). 
The sizes adopted in numerical calculations of five typical 
MPQDs are comparable to that in experiments^^. 


II. TIGHT-BINDING MODEL 

As shown in Fig. 1(a) , MLP has an irregular honeycomb 
structure with lattice constants <2=4.38A and h=?>3\A. There 
are four types of phosphorus atoms in a unit cell (Fig. 1(b)). 
The Hamiltonian for MPQDs can be written as^^ 
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where the summation runs over all the lattice sites of MPQDs, 
cj (Cj) is the creation (annihilation) operator of the electron at 
site i (j), and tij are the hopping energies. Five hopping links 
(see Fig. 1(a)) are needed to be taken into consideration^^. 
The related hopping integrals are ^ 1 =-1.220 eV, t2=3.655 eV, 
r3=-0.205 eV, r4=-0.105 eV, and r5=-0.055 eV. The band 
gap of MLP given by this TB model is 1.52 eV with the va¬ 
lence band maximum (VBM) and a conduction band mini¬ 
mum (CBM) located at -1.18 eV and 0.34 eV respectively^^. 
When we consider a magnetic field B applied perpendicularly 
to the plane of a MPQD, the transfer integral becomes tfj = 
^ where = | JI • A is the Peierls phase. In our 
calculation, we take the Landau gauge A = (0, Bx, 0) . The 
magnetic fiux S = ^ through a plaquette is in unit of Oq = |. 
The eigenvalues and eigenstates can be obtained from the sec¬ 
ular equation det\EI - H\=0, where . 

The electronic density of states (DOS) of a quantum dot 
(QD) are the sum of a series of Dirac S functions. Instead 
of plotting the level structure it is more useful to study the 
density of levels since it manifests more clearly the exact and 
nearly exact degeneracies of levels. Further, a regular varia¬ 
tion of the level density will indicate a shell structure. The 
energy DOS with a Gaussian broadening of a set of discrete 
levels 6i is given by 


D(€) = 




( 2 ) 


where F is the broadening factor, which is taken as 0.1 eV in 
our calculations without specification. 


TABLE I: Basic information of MPQDs including the index, number 
of phosphorus atoms Ntotai and the size of MPQDs studied. Na and 
are the number of hexagonal units along armchair and zigzag edge. 
RMPQDs are characterized by two indexs Na and N^, while other 
MPQDs are characterized by Na or N^. 



III. ELECTRONIC PROPERTIES OF MPQDS 

In this paper, five types of MPQDs will be explored in¬ 
cluding rectangular MPQDs (RMPQDs), hexangular armchair 


FIG. 2: (Color online) (a) Eigenenergies of the RMPQD with Na = 
Nz = 12. The X axis numbers the states with increasing energy, (b) 
DOS of the RMPQD. We use a Gaussian function with a broadening 
factor r=0.1eV to smooth the discontinuous energy spectra. The 
edge states are marked by red color. The horizontal green line is 
the Fermi level of this dot. 
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A. Rectangular MPQDs 

We begin with the electron structures in RMPQDs under 
zero magnetic field. The discrete levels and smoothing en¬ 
ergy DOS for a RMPQD with Na = = 12 are presented in 

Fig. 2. Edge states appear in the band gap of MLP (mid-in¬ 
gap states) below Fermi level {Ep) (red point in Fig. 2(a)). The 
corresponding DOS posses a peak within the MLP band gap 
(see Fig. 2(b)), which is below Fermi level {Ep) and marked 
by red color. The edge states play important roles in the elec¬ 
tronic and magneto properties of RMPQDs, which makes the 
RMPQDs different to MLP. The probability distributions con- 

(a) 

0.9 


0.8 


> 

- 3 ? 0.7 


>50.6 

0 


0.4 


0.3 

0 10 20 30 40 50 

Size 



(b) 



FIG. 3: (Color online) (a) The electronic density distributions of mid¬ 
in-gap states of RMPQD with Na - - 12. (b) Energy gap as a 

function of the size of RMPQDs. We change the dot size in three 
ways: (1) increase but fix A^=12; (2) increase Na but fix 12; 
(3) keep Na=Nz and increase them. Ecbm=0.34 eV is the conduction 
band minimum of monolayer phosphorene. 

tributed by all edge states in the RMPQD with = 12 

is shown in Fig. 3(a). We find that electrons in all edge states 
localize at the phosphorus atoms in the zigzag edge. Further , 
the number of mid-in-gap states is exactly equal to the number 
of phosphorus atoms in the zigzag edges in RMPQDs. This 
observation is important in tuning the electronic and magneto¬ 
optical properties of RMPQDs by changing the dot size. 

In Fig. 3(b), we examine the effect of the dot size (the pa¬ 
rameters Na and N^) on the energy gap for RMPQDs. We 
change the dot size in three ways: (1) increase but fix 
Na=l2; (2) increase Na but fix A^^=I2; (3) keep Na = N^ 
and increase them. In all three cases, the energy of the low¬ 
est unoccupied molecular orbital (LUMO) decreases with the 
dot size increasing, the energy of the highest occupied molec¬ 
ular orbital (HOMO) keeps constant (approximated by Ep). 
Therefore the gap between HOMO and LUMO does not ap¬ 
proach the energy gap of MLP (1.52 eV). When Na = N^, the 
energy of LUMO will decrease to the CBM of 2D phospho¬ 
rene. In the other two cases where either Na or N^ is fixed, 
the energy of LUMO decreases to a limit value higher than 
Ecbm of MLP, but approaches the CBM of corresponding rib¬ 
bon with a fixed width. When we increase the size by method 
(1), the LUMO and HOMO of the RMPQDs will approach 
the CBM and VBM of zigzag-edged monolayer phosphorene 
ribbons (ZMPR) with a width A^ = 12 respectively. The rea¬ 


son is that ZMPR poss topological protected edge states^^ . 
When we increase the size by method (2), only the LUMO of 
RMPQDs approaches the CBM of armchair-edged monolayer 
phosphorene ribbons (AMPR) with a width N = 12, while the 
HOMO of RMPQDs is higher than the VBM of AMPR with a 
width A = 12, because there are no edge states in the AMPR. 


B. Hexagonal and triangular MPQDs 

Now we explore the electronic properties of HAMPQDs, 
HZMPQDs, TAMPQDs, and TZMPQDs. For the four kinds 
of MPQDs, the numbers of hexagonal units (Na or N^) along 
all edges of their boundary are identical. 

In FIG. 4, we plot the DOS spectrum of electrons in hexag¬ 
onal and triangular MPQDs with different edge shapes. It is 
evident that mid-gap states appear near the Fermi energy Ep 
(see the red curves), regardless of the boundary shape or the 
edge types. The general existence of edge states in MPQDs 
is remarkable, which is different to that of graphene quantum 
dots (GQDs). In GQDs, edge states only appear in hexagonal 
zigzag edged GQDs^^’^^. 

The probability distributions of edge states in these MPQDs 
are shown in Fig. 5. For all edge states considered here elec¬ 
trons appear only at A-type and D-type (see Fig. 1(b)) ter¬ 
minated phosphorus atoms in the armchair direction, which 
is distinct to that in GQDs. In GQDs, the edge states are 
type-selected, i.e., the edges of a GQD are identical and the 
edge states appear or disappear in every edge of a GQD. How¬ 
ever, edge states in hexagonal and triangular MPQDs are not 
type-selected but direction-selected, that’s to say the electrons 
in edge states will appear at A-type and D-type terminated 
phosphorus atoms in the armchair direction , but disappear at 






FIG. 4: (Color online) Density of state of MPQDs. (a) HAMPQD, 
Na = 7. (b) HZMPQD, N, = 12. (c) TAMPQD, Na = 16. (d) 
TZMPQD, Nz = 20. The region where mid-gap states exist is marked 
by red color. 
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FIG. 5: (Color online) Electronic density distributions of edge states. 
(a)HAMPQD with Na = 7. (b) HZMPQD with N, = 12. (c) TAM- 
PQD with = 16. (d) TZMPQD with = 20. (e) Energy gap as a 
function of size of MPQDs. 


B-type and C-type terminated phosphorus atoms. This con¬ 
clusion is also suitable to RMPQDs. To illustrate this dif¬ 
ference, the easiest way is to find out the difference in the 
TB hamiltonian. In the GQDs, the hopping to three near¬ 
est neighbor carbon atoms are identical, which lead to the 
isotropic properties of GQDs. In MPQDs the hopping in¬ 
tegral to the nearest three phosphorus are ti (hopping along 
zigzag direction) and t 2 (hopping along armchair direction). 
The large difference between ti and t 2 makes MPQDs highly 
anisotropic. This is also the main reason why MPQDs have 
direction-selected edge state. 

The variation of the energy gap with the size of MPQDs 
is plotted in Fig. 5(e). The energy gap decreases as the size 
of the MPQD increases except for TAMPQDs. The energy 
of HOMO and LUMO of large size MPQDs are shown in 
TABLE. IT From TABLE. II, we can see that the HOMO 
of this four types of MPQD are from the edge states, while 
the LUMO are from the bulk state except for TAMPQDs. In 
TAMPQDs, both HOMO and LUMO are from the edge states. 
Therefore, for TAMPQDs the energy gap keeps constant as 
the dot size increases. 



HOMO(eV) 

LUMO(eV) 

Gap(eV) 

HAMPQD 

-0.11 

0.34 

0.45 

HZMPQD 

0 

0.34 

0.34 

TAMPQD 

-0.07 

0.10 

0.17 

TZMPQD 

0 

0.35 

0.35 


TABLE II: The energy of HOMO, LUMO and the gap between 
HOMO and LUMO when we increase Na (or N^) to 50. 


IV. MPQDS UNDER MAGNETIC EIELDS 

To reveal potential applications of MPQDs, it is necessary 
to study the magneto-energy spectrum and optical absorption. 
In this section we take the RMPQD with Na=Na=l2 as an 
example to examine the effect of magnetic fields on both the 



0/00 0/00 

PIG. 6: (Color online) (a) The DOS and energy spectrum of the RM¬ 
PQD with Na = Nz=l2 in the presence of magnetic field. We use 
a Gaussian function with broadening factor 0.01 eV here, (b) The 
magnetic energy spectra in relatively low magnetic field. The red 
lines correspond to the Landau levels of MLR 


energy spectrum and optical absorption. 

The calculated magneto-energy spectrum is shown in 
Fig. 6(a), which exhibits a clear Hofstardter-butterfiy feature. 
As the magnetic flux increases, the magnetic levels in the 
RMPQD approach the Landau levels of MLP^^, as shown in 
Fig. 6(b). A ’’fiat band” appears near the Fermi energy Ep, 
which is derived from the edge states at zero magnetic field. 
This is because the edge states have a negligible probability 
distribution with a distance smaller than or comparable to the 
magneto length (Ip = ^Jfi|eB). 

(a) (b) 
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FIG. 7: (Color online) (a) The optical-absorption spectrum a for 
RMPQD (Na = Nz = 12) in the absence of magnetic field. The 
black and red lines represent the absorption of light with polarization 
along the armchair and zigzag direction. We use a Gaussian function 
with broadening factor 0.005 eV to smoothen the discontinuous ab¬ 
sorption spectra, (b) The contour plot of the magneto-optical spectra 
of RMPQD (Na = A, = 12). 

The optical properties of MPQDs are promising for poten¬ 
tial applications in optic-electronic devices based on phospho- 
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rene. Therefore, we calculate the absorption spectra of RM- 
PQD with Na = Nz=\2, 

7ie^ 1 ^ |2 

a{fi(jj) = —r- > \s ' Pcv\ X d{Ec - Ey- fioj) 

m^socnojV ^ 

where n is the refractive index, c is the speed of light in vac¬ 
uum, So the permittivity of vacuum, mo is the free-electron 
mass, and ~s is the polarization vector of the incident light. 
The moment matrix is 

{n\ p \m) = imoin ^ ^ c* (?"' - r) {r'\H |r) 

r r' 

here {r'\H |r) is the hopping integral between r and r'. 

The optical-absorption spectrum for RMPQD with Na = 
Nz = 12 in the absence of magnetic field is shown in Fig. 7(a). 
The absorption of incident light with polarization along the 
armchair direction is about 20 times stronger than the absorp¬ 
tion of incident light with polarization along the zigzag direc¬ 
tion. The result is qualitatively the same to that of MLP^^. For 
incident light with polarization along the armchair direction, 
the absorption peak below 1.52 eV results from the transition 
from the edge states to the bulk states, which would be impor¬ 
tant in future applications of MPQDs. Note that the intensity 
of these absorption peaks are relatively weaker than those ab¬ 
sorptions due to transtion between bulk states. The main rea¬ 
son is that the phosphorus atoms on the MPQD boundary are 
much less than bulk phosphorus atoms. The effect of a mag¬ 
netic field on the optical spectrum of the RMPQD is shown 
in Fig. 7(b). In the contour magneto-optical graph, the strong 


absorption lines exhibit asymptotic behavior corresponding to 
the transitions between the conduction and valence band Lan¬ 
dau levels at a high magnetic field. 


V. SUMMARY 

We have theoretically investigated the electronic and 
magneto-optical properties of rectangular, hexangular and tri¬ 
angular MPQDs by means of the TB model. We find edge 
states detached from the bulk band gap in all kinds of MPQDs 
regardless of their shapes and edge configurations due to the 
anisotropic hoppings parameters in MLR For all edge states, 
the electron density is concentrated only in the armchair direc¬ 
tion of the dot boundary, which is distinct to that in graphene 
quantum dots. The magnetic levels of MPQDs exhibit a 
Hofstadter-butterfiy spectrum and approach the Landau lev¬ 
els of MLP as the magnetic field increases. A ’’fiat band” 
appears in the magneto-energy spectrum which is totally dif¬ 
ferent from that of MLP. The electronic properties and optical 
absorption of MPQDs can be tuned by changing the dot size, 
the types of boundary edges, and the external magnetic field. 
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